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Abstract 

This paper presents algorithms for computing the Grobner fan of an 
arbitrary polynomial ideal. The computation involves enumeration of all 
reduced Grobner bases of the ideal. Our algorithms are based on a uniform 
definition of the Grobner fan that applies to both homogeneous and non- 
homogeneous ideals and a proof that this object is a polyhedral complex. 
We show that the cells of a Grobner fan can easily be oriented acyclically 
and with a unique sink, allowing their enumeration by the memory-less re- 
verse search procedure. The significance of this follows from the fact that 
Grobner fans are not always normal fans of polyhedra in which case reverse 
search applies automatically. Computational results using our implemen- 
tation of these algorithms in the software package Gfan are included. 

1 Introduction 

The Grobner fan of an ideal I C k[xi, . . . , x n ] was defined by Mora and Rob- 
biano in 16 . It is a fan of polyhedral cones indexing initial ideals of I. The 
full-dimensional cones are in bijection with the distinct monomial initial ideals 
with respect to term orders or equivalently, the reduced Grobner bases of the 
ideal. In this paper we will describe algorithms for computing Grobner fans of 
arbitrary polynomial ideals allowing us to study their structure in detail. Our 
algorithms are implemented in the software package Gfan |12j . 

The computation of the Grobner fan of / in terms of reduced Grobner bases 
yields a universal Grobner basis of I, a set of polynomials which is a Grobner 
basis of I with respect to every term order. The Grobner fan also plays an 
important role in Grobner basis conversion [5] and the emerging field of tropical 
mathematics as it contains the tropical variety of I as a subfan [3] . Many of the 
well-known theoretical applications of Grobner bases rely on the existence of 
a Grobner basis or initial ideal for an ideal with prescribed properties such as 
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a particular complexity (a specified degree or squarefree-ness) or homological 
properties (Cohen-Macaulay, Gorenstein etc). The (even partial) computation 
of the Grobner fan makes such experimentations possible. No software package 
for the computation of Grobner fans is available at present. 

In the literature a distinction is often made between the case of I being 
homogeneous where the Grobner fan is a complete fan in W 1 , and the case of I 
being non- homogeneous, where the restricted Grobner fan with support equal 
to M>q is considered. In this paper we avoid this distinction by giving a new 
uniform definition of the Grobner fan and we prove that this actually defines a 
fan in the sense of a polyhedral complex — a proof that was left out in |16j but 
was proven for the special case of homogeneous ideals in JH]- See Section |2j 

In |18l Algorithm 3.2 and 3.6] and |161 Section 6] methods for computing 
the Grobner fan of a polynomial ideal were given. In this paper we shall study 
|18| Algorithm 3.6] in detail. This algorithm traverses the maximal cones of 
the Grobner fan. In the special case of toric ideals the traversal algorithm was 
already studied and implemented in . The traversal is graph-like — given a 
maximal cone we need to be able to find its facets and we need to be able to walk 
through a facet to the neighboring maximal cone. Algorithms for performing 
these local computations are discussed in Section |IJ These amount to solving 
linear programming problems and using the local basis change procedure due 
to jS] . We explain how to apply these methods to our case. 

The Grobner fan of a homogeneous ideal / is known to be the normal fan of a 
polytope, the state polytope of / ((T8J Theorem 2.5]). In this homogeneous case 
traversal of the maximal cones in the Grobner fan by walking through facets 
is equivalent to traversal of the edge graph of the state polytope. In the 
memory-less reverse search procedure for traversing the edge graph of a poly- 
tope was given. This procedure easily applies to Grobner fans of homogeneous 
ideals. However, the question is what happens if the ideal is not homogeneous. 
In JH] a non-regular Grobner fan was presented — a fan that is not the normal 
fan of any polyhedron. In light of this example it is not clear that the reverse 
search technique applies to Grobner fans in general. In Section |21 we prove that 
all Grobner fans have what we shall call the reverse search property, allowing 
them to be traversed efficiently. 

Grobner fans are often computed for ideals that possess a great deal of 
symmetry. In Section 14.41 we describe how to take advantage of symmetry in 
the computations. The methods used here are similar to those in Rambau's 
software package TOPCOM 17 for traversing the secondary fan of a point 
configuration up to symmetry. 

In Section |S] we discuss the complexity of our enumeration algorithm and in 
Section El we present several examples of Grobner fans computed using Gfan. 
This software package uses the GNU multi-precision library |1U| for exact arith- 
metics and Cddlib 7 j for solving linear programming problems. 
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2 The Grobner fan of a polynomial ideal 



Let R = k[xi, . . . , x n ] be the polynomial ring in n variables over a field k and 
let I C R be an ideal. The Grobner fan and the restricted Grobner fan of / 
are n-dimensional polyhedral fans defined in [TSj. We recall the definition of a 
fan in M n . A polyhedron in R n is a set of the form {x G K n : Ax < 6} where ^4 
is a matrix and b is a vector. Bounded polyhedra are called polytopes. If b = 
the set is a polyhedral cone. The dimension of a polyhedron is the dimension of 
the smallest affine subspace containing it. A face of a polyhedron P is either 
the empty set or a non-empty subset of P which is the set of maximizers of a 
linear form over P. We use the following notation for the face maximizing a 
form uj G W 1 : 

face^P) = {p G P : = maxg e p(w, <?)}. 

A face of P is called a /acei if its dimension is one smaller than the dimension 
of P. 

Definition 2.1 A collection C of polyhedra in R n is a polyhedral complex if: 

1. all non-empty faces of a polyhedron P G C are in C, and 

2. the intersection of any two polyhedra A, B G C is a face of A and a face 
of B. 

The support of C is the union of its members. A polyhedral complex is a fan if 
it only consists of cones. A fan is pure if all its maximal cones have the same 
dimension. 

A simple way to construct a fan is to take the normal fan of a polyhedron. 

Definition 2.2 Let P C W 1 be a polyhedron. For a face P of P we define its 
normal cone 

N P {F) := {w£M» : face w (P) = F} 

with the closure being taken in the usual topology. The normal fan of P is 
the fan consisting of the normal cones Np (P) as F runs through all non-empty 
faces of P. 

If the support of a fan is M n , the fan is said to be complete. It is clear that the 
normal fan of a polytope is complete. Not all fans arise as the normal fan of a 
polyhedron page 25]. 

For a G N n we use the notation x a := x" 1 . . . x^ n for a monomial in R. By 
a term order on R we mean a total ordering on all monomials in R such that: 

1. For all a G N n \{0} : 1< x a and 

2. for a, (3, 7 G N n : x° < ^ x a X T < x^. 
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By a term we mean a monomial together with its coefficient. Term orders are 
used for ordering terms, ignoring the coefficients. For a vector ui G M> and a 
term order -< we define the new term order as follows: 

x Q ^x^ (u,a) < (uj,(3) V ({u,a) = (u,P) A x a -< x' 3 ). 

Let -< be a term order. For a non-zero polynomial / G i? we define its initial 
term, in^(/), to be the unique maximal term of / with respect to -<. In the same 
way for ui G M. n we define the initial form, in w (/), to be the sum of all terms 
of / whose exponents maximize (uj, •}. The polynomial / is to -homogeneous if 
in w (/) = /. The uj-degree of a term cx a is (uj, a) and the uj-degree of a non-zero 
polynomial / is the common w-degree of the terms of in^(/). The initial ideals 
of an ideal I with respect to -< and uj are defined as 

in^(7) = (in x (/) : / G 7\{0}> and in,, (J) = (in^/) : / G I). 

Note that in^(J) is a monomial ideal while in^J) might not be. A monomial 
in i?\in^(7) (with coefficient 1) is called a standard monomial of in x (J). 

Although initial ideals are defined with respect to not necessarily positive 
vectors, Grobner bases are only defined with respect to true term orders: 

Definition 2.3 Let 7 C i? be an ideal and -< a term order on R. A generating 
set Q = {gi, . . . ,g m } for I is called a Grobner basis for I with respect to -< if 

in ^0O = (in^i), • • • M^(9m))- 

The Grobner basis Q is minimal if {in^(gi), . . . , m^(g m )} generates in^(I) min- 
imally. A minimal Grobner basis is reduced if the initial term of every g G Q has 
coefficient 1 and all other monomials in g are standard monomials of in^ (I) . 

We use the term marked Grobner basis for a Grobner basis where the ini- 
tial terms have been distinguished from the non-initial ones (they have been 
marked). For example, { x 2 + xy + y 2 } and {x 2 + xy + y 2 } are marked Grobner 
bases for the ideal (x 2 + xy + y 2 ) while {x 2 + xy + y 2 } is not since xy is not the 
initial term of x 2 + xy + y 2 with respect to any term order. 

For a term order -< and an ideal I, Buchberger's algorithm guarantees the 
existence of a unique marked reduced Grobner basis. We denote it by Q^(I). 
For two term orders -< and -<', if in^(7) = in^/(I) then Q^(I) = G^'(I). Con- 
versely, given a marked Grobner basis Q^(I), in^(I) can be easily read off. 

Given an ideal I, a natural equivalence relation on W l is induced by taking 
initial ideals: 

u ~ v in u (7) = m v (I). (1) 

We introduce the following notation for the closures of the equivalence classes: 

C^(I) = {u(£R n : m u (I) = in^(/)} and 

C V (I) = {u G R n : in u (I) = in„(/)}. 
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Remark 2.4 It is well known that for a fixed ideal I there are only finitely 
many sets C^(I) and they cover K> , see [TJI|- Secondly, every initial ideal 
in^(J) is of the form in^/) for some uj £ ^>o> see l^l Proposition 1.11]. 
Consequently, every C^(I) is of the form 0^(1). 

A third observation is that the equivalence classes are not convex in general 
since we allow the vectors to be anywhere in R n : 

Example 2.5 Let I = {x — l,y — 1). The ideal I has five initial ideals: (x — 
1, y — 1), (x, y), (x, y — 1), (x — 1, y) and (1). In particular, for u = (—1, 3) and 
v = (3, -1) we have in„(/) = in„(J) = (1) but ini (u+i;) (J) = (x,y). 

Proposition 2.6 Let -< be a term order and v € C^(I). For u € W 1 

in u (J) = in„(7) <==> V<? £ Q^(I), m u (g) = in v (g). 

This proposition is a little more general than Proposition 2.3 in ^H] as it allows 
the vectors u and v to have negative components. A proof is given in the next 
section. For fixed -< and v as in Proposition l2.61 we get that C V (I), the closure of 
the equivalence class of v, is a polyhedral cone since each g € 0^,(1) introduces 
the equation in u (g) = in v (g) which is equivalent to having u satisfy a set of 
linear equations and strict linear inequalities, see Example 12.71 The closure is 
obtained by making the strict inequalities non-strict. Under the assumptions 
of Proposition 12.61 we may write this in the following way: 

ueC v (I) Vg £ g^(I), m v (in u (g)) = in v (g). (2) 

As we saw in Example 12.51 not all equivalence classes are convex. However, 
for an arbitrary v, C V (I) is a convex polyhedral cone if it contains a strictly 
positive vector. In this case, there must exist a vector p £ in the interior 
of C V (I). Then in p (J) = m v (I) and, by Lemma f2.151 p £ C^ p (I) for any -<. 
Hence the equivalence class of v is of the form required in Proposition 12.61 

Example 2.7 Let / = (x + y + z, x 3 z + x + y 2 ) C Q[x, y, z] and let -< be the 
lexicographic term order with x ~< y -< z. Then G~>(I) = {y 2 + x — x^y — x 4 , z + 
y + x}. Ev = (1,4,5) then in^I) = in^(J) = {y 2 ,z) a,ndC v (I) = C^(I). By 
Proposition 12.61 in u (J) = in v (I) if and only if the following two equations are 
satisfied: 

in u (z + y + x) = z (•<=>• u z > max{n x , u y }), and 

in u (y 2 + x - x 3 y - x 4 ) = y 2 2u y > m&x{u x , 3u x + u y ,4u x }). 

Introducing non-strict inequalities we obtain a description of C^(I). This 
cone is simplicial and has the cones CV 0) o,i)(-0) Ql,3,3)C0 an d C(-2,-l,-l) (-^ ) 
as extreme rays and Cn )3) 4)(J), C/ 2,-1,0) (-0 an d C(-i,2,2) (-0 as facets. Since 
(—2, —1,0) is in C^(I) a description of vectors u in C(_ 2 ,-i,o)(-0 i s given by: 

m (-2,-i,o)( in «( 2 + U + z)) = z -u z > max{n a; ,u y }), and 

m (-2-i,o)( in M(y 2 + x - a:3 2/- :r4 )) = y 2 +^ 2u y = u x > max{3u x + u y ,4u x }). 
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Figure 1: The Grobner fan of the ideal in Example 12 . 71 has 7 three-dimensional, 
14 two-dimensional and 8 one-dimensional cones. The intersections of the two- 
dimensional cones with the hyperplane x + y + z = 1 are drawn as lines. The 
dotted part of the figure shows the combinatorial structure outside the hyper- 
plane. The gray triangle indicates the positive orthant. 

Definition 2.8 The Grobner fan of an ideal / C R is the set of the closures 
of all equivalence classes intersecting the positive orthant together with their 
proper faces. 

This is a variation of the definitions appearing in the literature. The advantage 
of this variant is that it gives well-defined and nice fans in the homogeneous and 
non-homogeneous case simultaneously. By nice we mean that all cones in this 
fan are closures of equivalence classes. It is not clear a priori that the Grobner 
fan is a polyhedral complex. A proof is given in the next section (Theorem 
12.19(1 . The support of the Grobner fan of / is called the Grobner region of /. 
Recall that the common refinement of two fans J-\ and Ti in W 1 is defined as 

T\ A T 2 = {C\ l~l C2}(c 1 ,C 2 )€3 r 1 xF 2 - 

The common refinement of two fans is a fan. We define the restricted Grobner 
fan of an ideal to be the common refinement of the Grobner fan and the faces 
of the non-negative orthant. The support of the restricted Grobner fan is R>n. 
The Newton polytope of a polynomial is the convex hull of its exponent vectors. 

Example 2.9 The Grobner fan of the principal ideal (x 4 + x 4 y — x^y + x 2 y 2 + 
y) consists of one O-dimensional cone, three 1-dimensional cones and two 2- 
dimensional cones, see Figured The same is true for the restricted Grobner fan. 
Notice, however, that in the restricted Grobner fan one of the 1-dimensional 
cones and one of the 2-dimensional cones are not equivalence classes of the 
equivalence relation ©. 
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Figure 2: The Grobner fan of the ideal in Example 12.91 is shown on the left. 
The restricted Grobner fan is on the right. In the middle the Newton polytope 
of the generator is drawn with the shape of its normal fan indicated. 

2.1 Proof: The Grobner fan is a fan 

In this section we prove that the Grobner fan is a fan i.e., that it is a polyhedral 
complex consisting of cones. Recall, in general the Grobner fan is not complete 
and its support is larger than M> . In |16| there is no proof that the Grobner 
fan is a fan in the sense of a polyhedral complex. A proof that the Grobner fan 
is a polyhedral complex under the assumption that the ideal is homogeneous is 
given in ^Sj. We present a complete proof for the general case. Many of the 
results we need in the proof are generalizations of known results needed in the 
proof that the Grobner fan of a homogeneous ideal is a polyhedral complex |18j . 
However, we do not rely on these references for the sake of being self-contained. 

We fix the ideal I C R in the following theorems. The most important step 
is the proof of Proposition 12.61 which tells us that the closure of an equivalence 
class is a polyhedral cone. Then we prove that the relative interior of any face 
in the Grobner fan is an equivalence class (Proposition I2.l5)l and, finally, that 
the intersection of two cones in the fan is a face of both (Proposition 12. . 

To prove Proposition 12.61 we start by proving a similar statement for the 
equivalence classes arising from initial ideals with respect to term orders. 

Lemma 2.10 Let -< be a term order. For v £ W 1 , 

ULu(J) = in^(I) <=4> V<7 £ Q^(I), in v (g) = 'm^(g). 

Proof. =4>: Let g £ G^(I). Since Q^(I) is reduced, only one term from g, in^(^), 
can be in in^(/) = in„(J). The initial ideal in„(I) is a monomial ideal, implying 
that all terms of an element in the ideal must be in the ideal too. Hence, the 
initial form in v (g) £ in„(J) has to be equal to in^(<7). 

-<=: We must show that in„(I) = 'm^(I) where in_<(i") = (in^(<7)) ge g_ < (/). The 
"5" inclusion is clear since in^(g) = in v (g) £ in„(J) for all g £ Q^I). 

To prove the "C" inclusion, since in„(I) = (in„(/), / £ I), it suffices to 
show that irLy(/) £ in^(J) for all / £ /. Pick / £ I and reduce it to zero using 
the division algorithm (e.g. 6, Chapter 2]) with Q^{I) and -<. We may write 

/ = rnig^ H h m r g ir (3) 
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where rrij is a monomial and Q{. is an element from Q^{I). The division algo- 
rithm guarantees that in^(/) > m 3 -in_< (g^. ) with respect to -< since monomials 
are substituted with monomials less than the original ones with respect to -< in 
the division process. Exactly the same thing is true for u-degrees since v and 
-< agree on Q^(I). Thereby, any monomial on the right hand side in © has 
■u-degree less than or equal to the f-degree of the left hand side. Consequently, 

m f(/) = ^2mjm v (g i:j ) 
jeJ 

with j running through a subset such that mjin v (gi ) has the same u-degree as 
m v (f). Since in v (g) G in^(7), the initial form m„(/) G in._<(J). □ 

By Lemma 12.101 the equivalence class of in^(J) is open. Since in_;(I) is of 
the form hxy(I) for some v (see Remark 12.4(1 . the equivalence class of in^(I) 
is also non-empty and hence full-dimensional. Thus we have proved that the 
equivalence class of a term order is a full dimensional open polyhedral cone. 

Corollary 2.11 Let -< be a term order and v G W 1 . Then 

v G C y (J) O V# G g^(I) : m^m^g)) = in^(g). 

Proof. Lemma 12.101 tells us that v lies in the interior of (I) if and only if 
in v (<7) = in^(g) for all g G G^(I). Relaxing the resulting strict inequalities to 
non-strict inequalities we get a description of C^(I). This relaxation is exactly 
the one given by in^(in„(<?)) = in^(^) for all g in Q^(I). □ 

Lemma 2.12 A polynomial f G in„(J) can be written in the form f = in v (cj) 
where Cj G I and all summands in the sum have different v -degrees. 

Proof. The initial ideal in v (I) is generated by v- homogeneous polynomials, im- 
plying that all ^-homogeneous components of / are in in„(/). Let h be a maxi- 
mal ^-homogeneous component of /. We need to show that h is the initial form 
of an element in / with respect to v. We may write h as in„(ai) + • • • + in„(a s ) 
for some polynomials a\, . . . ,a s in I. Since h is v- homogeneous we can rewrite 
h as the sum YljeJ in„ (a^ ) of forms having the same f-degree as h. We pull out 
the initial form and get h = ^v(Ylj^j a j)- ^ 

Lemma 2.13 Let -< be a term order. If v G C^(I) then in-<(iau(J)) = in^(J). 

Proof. Let g G Q^(I). Since v G C^(I), by Corollarv l2.111 in^(g) = in^(in„(p)) 
and hence in^(J) = (in x (flO)<ye^(/) ^ "^(in^I)). 

We now prove that in^ (in„ (I) ) C in^ (I) . Notice that in_< (in,, (I) ) is gener- 
ated by initial terms of elements / G in v (/)\{0} with respect to -<. Suppose 
/ G in„(J)\{0}. It suffices to show that in^(/) G in^(7). Using Lemma 12.121 
we may write / = ^i=i m «( c «) where ci,...,c s G / and m v (ci), . . . ,in v (c s ) 
are u-homogeneous each with distinct degree, so that no cancellations occur. 
Consequently 'm^(f) equals in^(in„(cj)) for some j. We wish to prove that 
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in^(in„(cj)) € in^(J). We use the division algorithm with G^{I) and -< to 
rewrite c,- 

cj = miSd H h m r 5f ir 

where mi, . . . , m r are monomials and g^,. . . ,gi r belong to Q^{I). Let M be 
the v-degree of Cj. In the division algorithm we sequentially reduce Cj to zero. 
In each step, the u-degree of Cj will decrease or stay the same since we subtract 
the product of a monomial and an element from where the f-degree of the 

product already appeared in Cj by Corollary 12.111 Equivalently, the product of 
the monomial and the element from Q^{I) are "added" to the right hand side 
of the equation. We are done when c, = and or equivalently, the original Cj is 
written as the above sum with every term having u-degree less or equal to M. 
Consequently, we have 

in v (cj) = ^2 m «("V5v) 
j'eJi 

for a suitable J'. The division algorithm guarantees that the exponent vec- 
tors of in^fmigjj), . . . , m^(m r gi r ) are distinct. Since v G C^(J), they equal 
in_<(ra (mi<7i 1 )), . . . ,'m^(m v (m r gi r )). The maximal one of these with respect 
to -< cannot cancel in the sum. Hence in^(m v (cj)) = ia^(Tnjig% , ) for some f 
which implies that in^(in„(cj)) G in^(I) as needed. □ 

An easy corollary is a method for computing Grobner bases for initial ideals. 

Corollary 2.14 Let -< be a term order. If v G C^(I) then 

g^(m v (I)) = {m v (g)} geg ^ (I) . 

Proof. By C orollary |2JJJ (m^(m v (g))) geg ^ {1) = (in^(g)) geg ^ {1) = in^(J). By 
Lemma l2.13l in^(I) equals in^fin^f/)). Thus in^(in„(/)) = (im^in^)))^^/). 
This proves that {^v(g)} g ^g^(i) is a Grobner basis of in^-T) with respect to -<. 
It is reduced since Q^{I) is minimal and reduced. □ 

We are now able to give a proof for Proposition 12.61 which claimed that given 
v G C^(I) and u G W 1 , in u (7) = in„(J) <=^ V# G 9^(1), in u (g) = in v (g). 

Proof. <=: Since in u (g) = ixi v (g) for all g G Q^(I), we get that in^(in u (<?)) = 
in^(in„(5f)) for all g G Q^(I). Since v G C^(I), by Corollary 12.111 'va^(g) = 
m^(m v (g)) for all g G G-<{I) and hence in^(g) = in^(in u (g)) for all g G G-<{I) 
and u G C-^(J) by Corollary 12.111 The Grobner basis Q^(in u (I)) is then 
{^u(g)}geQ^(I) by Corollary 12.141 We get the same Grobner basis for in v (I). 
Hence, in u (I) = in„(I). 

=>: Let g G Q^{I). We need to show that \a u {g) = 1^(5). Since the basis is 
reduced, only one term of g, namely in^(g), is in in^(J). We start by proving 
that the term m^(g) is a term in in v (g) and a term in m u (g). For in«(<7) we apply 
Corollary 12.111 which says in^(g) = 'm^(in v (g)) . For in u (g) we apply Lemma 
12.131 and get in^(in u (p)) G in^(in u (/)) = in^(in w (J)) = in^(J). Only one term 
of g is in in^(J), so m^(in u (g)) = in^(g). If the difference in u (g) — in v (g), 
belonging to in u (7) = in^(I), is non-zero we immediately reach a contradiction 
since the difference contains no terms from in^(I) = in^(in„(/)). □ 
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We have now proved that every equivalence class of a vector v in a (/) is a 
relatively open convex polyhedral cone. By the argument following Proposition 
12.61 in the previous section all sets in the Grobner fan are in fact cones. We 
now argue that the relative interior of every cone in the Grobner fan is an 
equivalence class. 

Lemma 2.15 Let -< be a term order. If v £ M> then v £ C^ V (I). 

Proof. This follows from Corollary 12.111 since in^ v (in v {g)) = in^ v (g) for all 

Proposition 2.16 The relative interior of a cone in the Grobner fan is an 
equivalence class (with respect to u ~ u' -v4> in u (I) = m u >(I)). 

Proof. By definition every cone in the fan is the face of the closure of an equiv- 
alence class for a positive vector v £ K™ . Let -<' be an arbitrary term order 
and define -< as -<' v . According to Lemma 12 . 1 51 the vector v belongs to C^(I). 
Notice that by ©, C V {I) C C X (J) since for all u £ C V (I) and g £ G~<(I), the 
cond ition in^(in u (sr)) = m^(m v (m u (g))) = m^('m v (g)) = in^(g) of Corollary 
12.111 is satisfied. By (J2J) the closed set C V (I) is cut out by some equations and 
non-strict inequalities. The relative interior of any face of C V (I) can be formed 
from this inequality system by changing a subset of the inequalities to strict 
inequalities and the remaining ones to equations. So let u be a vector in the 
relative interior of some face of C v (I) . The vector u is in C v (I) C C_< (I) . We 
may use Proposition 12.61 to conclude that a vector u' £ W 1 is equivalent to u if 
and only if it satisfies the inequality system mentioned above — that is, if and 
only if it is in the relative interior of the face. □ 

It remains to be shown that the intersection of two cones in the Grobner 
fan is a face of both cones (Proposition 12. l8]l . We need a few observations. 

Corollary 2.17 Let C be a cone in the Grobner fan. Ifv £ C then for u £ W 1 , 

in u (J) = \n v {I) => u £ C. 

Proof. The vector v is in the relative interior of some face of C. This face is 
also in the Grobner fan. By Proposition 12. 161 u is in the relative interior of the 
same face and, consequently, also in C. □ 

By Remark 12.41 there are only finitely many initial ideals given by term 
orders and, consequently, only finitely many reduced Grobner bases of /. It 
follows that there can only be finitely many equivalence classes of the type 
described in Proposition 12.61 and Proposition 12.161 

Proposition 2.18 Let C\ and C2 be two cones in the Grobner fan of I. Then 
the intersection C\ H C2 is a face of C\ . 
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Proof. The intersection C\ n C2 is a cone. By Corollary 12.171 Ci and C2 are 
unions of equivalence classes. Further, if v G Ci n C2, then again by Corollary 
12.171 the entire equivalence class of v is both in C\ and in C2 and hence in 
C\ fl C2. Hence Ci n C2 is a union of equivalence classes. 

Let ii be a vector in such an equivalence class E contained in C\ fl C2 . Then 
■u is in the relative interior of one of the faces of C\ which is a cone in the 
Grobner fan. By Proposition 12.161 the set of vectors in the relative interior of 
this face is exactly E. Hence every such equivalence class is the relative interior 
of a face of C\ and its closure is the face. 

Look at the R-span of each equivalence class contained in C\ n C2. These 
spans must be different for every face of C\. We claim that there can be only 
one maximal dimensional cone/span. If there were two cones then their convex 
hull would be in C\ n C2 and have dimension at least one higher and thus 
cannot be covered by the finitely many lower dimensional equivalence classes 
— a contradiction. 

Let E be the maximal dimensional equivalence class contained in C\ fl Ci- 
We will argue that E = CiflC^. The inclusion E C C\^C2 is already clear since 
C\ fl C2 is closed. To prove the other inclusion suppose u> G C\ fl Cz\E, Then 
conv(E,uj)\E is contained in C\ IIC2 and has dimension at least the dimension 
of E. This is a contradiction since conv(-E, u)\E cannot be covered by finitely 
many lower dimensional equivalence classes. This completes the proof. □ 

Theorem 2.19 The Grobner fan is a polyhedral complex of cones and hence a 
fan. 

Proof. We already argued using Proposition 12.61 and Lemma 12.151 that the 
Grobner fan consists of polyhedral cones. The first condition for being a poly- 
hedral complex is satisfied by definition. The second condition is Proposition 

EH □ 

3 Reverse search property 

By the graph of a pure full-dimensional fan we mean the set of maximal cones 
with two cones being connected if they share a common facet. In this section 
we will prove that the reverse search technique can be used for traversing 
the graph of a Grobner fan. This follows from the main theorem, Theorem 
.61 which says that the graph of a Grobner fan can be oriented easily without 
cycles and with a unique sink. In Definition 13.41 we define what we mean by 
this. 

We start by explaining how a graph with this special kind of orientation 
can be traversed by reverse search. The idea is to define a spanning tree of 
the graph which can be easily traversed. The following is a simple proposition 
which we shall not prove. 

Proposition 3.1 Let G = (V,E) be an oriented graph without cycles and with 
a unique sink s. If for every vertex v G some outgoing search edge 

e v = (v, •) is chosen then the set of chosen edges is a spanning tree for G. 
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The spanning tree in Proposition 13.11 is referred to as the search tree. The 
proposition implies that the graph is connected. 

Notice that we can find the sink by starting at any vertex and walking 
along a unique path of search edges until we get stuck, in which case we are at 
the sink. Consequently, the sink is the root of the oriented spanning tree. A 
corollary to the proposition is the reverse search algorithm for traversing G: 

Algorithm 3.2 Let G = (V, E) be the oriented graph of Proposition \3.1\ and 
suppose the choice of a search edge e v for each vertex v ^ s has been made. 
Calling the following recursive procedure with v = s will output all vertices in 
G. 

Output .subtree (v) 

Input: A vertex v in the graph G. 

Output: The set of vertices in the subtree with root v. 

{ 

Output v; 

Compute the edges of form (-,v) £ E; 
For every oriented edge (u, v) £ E 

If ( e u = ( u i v )) Output_subtree^j; 

} 

This algorithm does not have to store a set of "active" vertices as is usually 
needed in depth- and breadth-first traversals. It is even possible to formulate 
the algorithm completely without recursion avoiding the need for a recursion 
stack. In that sense the algorithm is memory-less. 

We give an example of how the edge graph of a polytope or, equivalently, 
the graph of its normal fan can be oriented. 

Example 3.3 Let P C W 1 be a polytope whose vertices have positive integer 
coordinates and let -< be a term order on R. The following is an orientation of 
the edge graph of P without cycles and with a unique sink: An edge (p, q) is 
oriented from p to q if and only if x p -< x q . 

This defines an orientation of the graph of the normal fan of a polytope 
for any term order. We would like to mimic this orientation for any pure full- 
dimensional fan in M n . For simplicity we shall restrict ourselves to fans whose 
(n — l)-dimensional cones allow rational normals. In view of Propositions 12.61 
and 12.161 this is no restriction for Grobner fans. 

Definition 3.4 A pure full-dimensional fan in M. n is said to have the reverse 
search property if for any term order -< the following is an acyclic orientation 
of its graph with a unique sink: If {C\,C2) is an edge then C\ and Ci are 
n-dimensional cones with a common facet F. Let p,q £ N n such that q — p ^ 
is a normal for F with all points in C\\F having negative inner product with 
q — p and all points in C%\F having positive inner product with q— p. We orient 
the edge in direction from C\ to C2 if and only if x p -< x 9 . 

Note that the orientation of an edge in Definition 13.41 does not depend on the 
particular choice of p and q. Note also that for normal fans of polytopes this 
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Figure 3: A fan not having the reverse search property, see Example 13,51 



orientation agrees with the orientation of the edge graphs of the polytopes in 
Example 13.31 Not every fan has the reverse search property: 

Example 3.5 Figure 01 shows a fan with support M> intersected with the 
standard simplex. The intersection is the non-dotted part of the figure. For 
every shared 2-dimensional facet the orientation of its edge with respect to a 
term order of form -<ni t i-\ is indicated by an arrow. The graph has a cycle. The 
reason is that the vector (1,1,1) is in the interior of the cone over the dotted 
triangle and therefore induces the shown orientation with any tie-breaking. 

Example 13.31 on the other hand shows that any normal fan of a polytope 
has the reverse search property. If I is a homogeneous ideal the Grobner fan of 
/ is known to be the normal fan of the state polytope of /, see ^5] for a proof. 
(We should mention that in |15] it was proven that this is only true if we use 
the state polytope definition in |18j and not true with the original definition in 
[2].) As a consequence the Grobner fan will have the reverse search property. 
The reverse search orientation of a fan with respect to any term order can be 
carried out on any fan covering M> and being the normal fan of a polyhedron. 
Since the restricted Grobner fan of any 0-dimensional or principal ideal satisfies 
these conditions it is clear that these fans have the reverse search property. 

It is shown in ^3] that this line of reasoning cannot be applied to Grobner 
fans in general. In particular, an ideal is presented whose restricted Grobner 
fan is not the normal fan of a polyhedron. For this reason we need a non-trivial 
argument to prove the following theorem: 

Theorem 3.6 The Grobner fan of any ideal I C R has the reverse search 
property. 

The proof is given in the next section. In Section 0] we will argue that all parts 
of Algorithm 13.21 (finding adjacent edges, finding adjacent vertices and finding 
search edges) can be implemented efficiently for Grobner fans. 
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3.1 Proof: The Grobner fan has the reverse search property 

In this section we prove Theorem Yd .61 We start by recalling how the polynomial 
ring can be graded by semigroups. This leads to a more general notion of 
homogeneous ideals. 

Definition 3.7 By a grading on R = k[x\, . . . ,x n ] we mean a pair (A, A) 
consisting of an abelian semigroup A and a semigroup homomorphism: 

A-.n n ^a 

such that A -1 (a) is finite for all a 6 A. The A- degree of a term cx b is A(b). A 
polynomial is A-homogeneous if all its terms have the same ^-degree. An ideal 
is A-homogeneous if it is generated by a set of .4-homogeneous polynomials. 

For a grading (A, A) on R we get the direct sum of /c-vector spaces 

where R a denotes the fc-subspace of R consisting of ^-homogeneous polynomials 
of degree a. Any reduced Grobner basis of an „4-homogeneous ideal I consists 
of ^4-homogeneous polynomials. In particular, by generalizing the argument of 
Lemma 12.121 we get the direct sum 

where I a denotes the A;-subspace of I consisting of ^-homogeneous polynomials 
of degree a. The ^4-homogeneous part I a is a fc-subspace of R a . We define the 
A- graded Hilbert function: 

H IA :A -» N (4) 
a i-> dim k (R a /I a ) (5) 



Remark 3.8 For a monomial ideal / the standard monomials of degree a form 
a basis for R a /I a . Hence Hj^(a) counts the number of standard monomials of 
degree a. 

In general, as the following well-known proposition shows, the Hilbert function 
can be found by looking at a monomial initial ideal: 

Proposition 3.9 Let I be an A-homogeneous ideal and -< a term order then 

Hi, a = #in^(J),.4- 

Proof. The linear map taking a polynomial to its unique normal form by the 
division algorithm on Q^{I) induces an isomorphism of fc-vector spaces 

Ra/Ia -> Ra/bu(I) a . 

□ 
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Consider a shared facet of the cones C\ and C 2 in the Grobner fan with a 
relative interior point v. The "edge ideal" in v (I) is homogeneous with respect 
to any vector in the relative interior of the facet and consequently also homoge- 
neous with respect to any vector in the span of the facet. Since C\ and C 2 both 
contain positive vectors, so does span K (C„(/)). Recall that C V (I) is the closure 
of the equivalence class of v. Pick a basis u\, . . . ,u n -i € N n for span R (C„(/)) 
with u\ being a positive vector. The vectors induce a grading A v : N n — > N n_1 
on R by 

A v (b) = ((ui,b), . . . , 6)) 
for 6 € N n . The initial ideal in„(J) is ^-homogeneous. 

Lemma 3.10 Let -< be a term order, I an ideal, (Ci,C 2 ) a directed edge with 
respect to the orientation in Definition ^. J\ and M\ and M 2 the initial ideals of 
C\ and C 2 respectively. Let v be a relative interior point in the shared facet. 
Then in^(in„(J)) = M 2 . 

Proof. Choose a positive interior point u 2 of C 2 . We claim that the following 
identities hold: 

M 2 = m W2 (L) = in^^in^I)) = in^ 2 (J) = in^ W2 (in c (J)) = m^(m v (L)). 

The first one holds by the choice of UJ2- The second one is clear since in^ (/) is 
a monomial ideal. The third one holds by Lemma 12.131 and Lemma 12.151 By 
Lemma \2. 131 the fourth equality holds since v € (I) = C a;2 (I). To prove 
the last equality we look at the reduced Grobner basis Q^(m v (L)). If we can 
show that in^ 2 (g) = m^(g) for all elements g £ Q^('m v (I)) then we know 
that ^_<(in^(/)) is also a Grobner basis with respect to and the generators 
for the initial ideal in^ (in„(I)) are exactly the same as those for in^(in„(/)). 
This would complete the proof. 

The reduced Grobner basis Q^(in v (I)) is ^-homogeneous. For an element 
g this implies that the difference between two of its exponent vectors must be 
perpendicular to the shared facet. By Definition 13 .41 there exists a normal q — p 
of the facet with x p -< x 9 and (u>2iQ — p) > 0. Since -< and agree on 
one normal vector they must agree on all exponent differences of elements in 
g^(m v (L)). □ 

Notice that by Proposition 13. 91 anv initial ideal in^(in„(/)) of in„(J) has the 
same ^-graded Hilbert function as in v (I). 

By a flip we mean a move from one vertex in the graph to a neighbor. 
For a degree a 6 N n_1 we call A~ 1 (a) the fiber over a. The ^-graded Hilbert 
function of an initial ideal in^in^I)) counts the number of standard monomials 
inside each fiber. A flip preserves the Hilbert function. We may think of this as 
monomials in the monomial initial ideal moving around in the fiber. We wish 
to keep track of how the monomials move when we walk in the oriented graph. 
We define exactly what we mean by "moving around" : 
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Definition 3.11 Let -<, Mi, M2, ui, ■ ■ ■ , u n -\ and v be as above with 
in^(in„(/)) = M2. Let N\ and N2 be the monomials in Mi and M2 respec- 
tively. We define the bijection 4>^M\hh '■ N\ — > N 2 in the following way: For 
a monomial x b 6 TVx look at the monomials £?i Q Ni and i?2 C N2 with the 
same .4-degree as x 6 . Since taking initial ideals preserves the ^4-graded Hilbert 
function, |jBi| = \B 2 \. Sort Bi and B2 with respect to -<. The bijection 4>^m 1 M2 
is now defined by taking the first element of Bi to the first element of B2, the 
second element of B\ to the second element of B2 and so on. 

The following lemma is from |14[ Lemma 4.1]: 

Lemma 3.12 Let <i and < 2 be two term orders. If , . . . , fg is a vector space 
basis for I a such that in< 1 (/-j 1 ), . . . , in< 1 (/g) is a basis for in< 1 (J) a , then there 
exists a basis ff, ■ ■ ■ , fg f or la such that in< 2 (/^), . . . ,in< 2 (/^) is a basis for 
in< 2 (I) a and 

n< 3 {ft)<in<Mi) 

in< 2 (/ s 2 ) <iin< x {fl). 
Corollary 3.13 Let the setting be as in Definition IS. Ill If x fe G Mi then 

Proof. Let a be the .4-degree of x b . We apply Lemma 13.121 with / in the 
lemma being in„(J). Let <i be -< and <2 be the refinement of the pre- 
order induced by ui with the reversed order of -<. By the orientation of the 
graph Mi = in< 2 (in„(/)) and M2 = in< 1 (in„(/)). By multiplying elements of 
Q^(in v (I)) by monomials we can construct a A:-basis f^, . . . , fj of in v (/) a with 
m.< il {f\), . . . , in< 1 (/g) being a basis of (M2) a . By the lemma there is a basis 
in< 2 (/j), . . . ,in< 2 (/ 2 ) of (Mi) a . Sort the list of inequalities in the lemma with 
in< 2 (/?) decreasing w.r.t. -< (<i). The right hand side can now be sorted with 
respect to the same order without violating the inequalities. To see this use the 
bubble sort algorithm — when two adjacent inequalities are swapped . . . 

in< 2 (/2) <i in^/i) in< 2 (/2) <, m^+i) 

Vi Ai i-» Vi Vi 

in < 2 (/m) <i i n <i(//+i) in < 2 (/m) <i '^<i(f!) 

. . . the relations on the right hand side of the arrow hold by transitivity of <i. 

After sorting, x fe appears somewhere on the left and 0-<AfiAf2( x ) on the 
right in the same inequality. This completes the proof. □ 

Proof of Theorem \S.(\ Suppose C\, C 2 , ■ ■ ■ , C m was a path in the oriented graph 
with Ci = C m . Let Mi, . . . , M m denote the initial ideals and Ni, . . . , N m their 
monomials. We will prove that the bijection <f> := (p^M m - 1 M m o ■ • • o 4>^m-i_m 2 ls 
the identity on M\. Suppose it is not the identity and let x 6 be the smallest 
element in Mi with respect to -< that is not fixed by <f>. By Corollary 13.131 x b 
is the image of a smaller element in Mi with respect to -<. But this element 
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is fixed by the minimality of x b — a contradiction. The composition being the 
identity implies by Corollary 13 1 131 that 4>^M t M t+1 is the identity for all i. Hence 
Mi = M i+ i, contradicting that M 1; M 2 , . . . , M m is a path. 

We claim that C^(I) is the unique sink. If v is in the relative interior of a 
facet of C^(I) then by Lemma 12. 131 in^fim.f ID = in^(I). By Lemma 13. 1UI this 
means that all edges connected to C^(I) are ingoing. Hence C-^(J) is a sink. 

To prove uniqueness let C^'(I) be some sink in the oriented graph. By 
|16j -< has a matrix representation (to, . . . , T n _i) € M nxn such that r e := tq + 
sti + • • • + e™ _1 r„_i € int C^(I) for e > sufficiently small. Furthermore, for 
any / € it!, in Te (/) = in^(/) for e > sufficiently small. If C^'(J) is a sink 
then according to Definition 13.41 there exists a complete list of inner normals 
qi — pi, . . . , q r — p r of C^i(I) fl M> such that in_<(x 9i — x p,; ) = x 9 \ Since r e 
and -< pick out the same initial forms on a finite set of polynomials for e > 
sufficiently small we see that (r £ ,qi) > (r £ ,pi) or, equivalently, t £ S int C_</(J) 
for e > sufficiently small. We conclude that C^/(J) = C_<(/). □ 

4 Implementation issues 

We can find a single Grobner cone by applying Buchberger's algorithm and 
Corollary 12. 1 II for some term order. Since the graph of the Grobner fan of I is 
connected we may choose any graph traversal algorithm for computing the full 
dimensional Grobner cones. To do the local computations we need to be able 
to find the edges (connecting facets) of a full dimensional cone and we need to 
be able to find the neighbor along an edge. We will see how to do this in the 
following sections. 

Throughout the graph enumeration process we will represent the Grobner 
cones by their marked reduced Grobner bases, rather than by their defining 
inequalities, their term orders etc.. This choice is justified by the following 
known theorem which we shall not prove: 

Theorem 4.1 Let I C R = k[x\, . . . ,x n ] be an ideal. The marked reduced 
Grobner bases of I, the monomial initial ideals of I (w.r.t. a positive vector) 
and the full- dimensional Grobner cones are in bijection. 

An important issue when implementing the algorithms is to identify shared 
facets. We say that a facet is flippable if its relative interior contains a positive 
vector. The flippable facets in a Grobner fan are always shared. With the right 
definition of search edges the search tree will only consist of flippable facets. 

At the end of the section we will see how the search edge computation in 
the reverse search algorithm can be implemented and we will explain how to 
take advantage of symmetry in a Grobner fan traversal. 

4.1 Finding facets 

Suppose that we know a marked reduced Grobner basis G~<(I) with respect to 
some unknown term order -<. Proposition 12.61 for Corollary 12. tells us how 
to read off the defining inequality system for C^(J). 
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Since C^(I) is full-dimensional the system contains no equations but only 
inequalities. Some of these inequalities are equivalent in the sense that they are 
multiples of each other. Taking just one inequality from each equivalence class 
the problem is now to find irredundant facet normals of a cone — or equivalently 
to find the extreme rays of the dual cone. Checking if a ray is extreme can be 
done by linear programming. 

Not all of the remaining inequalities are guaranteed to define flippable facets. 
One way to ensure that we only get flippable facets is by adding the constraints 
ej • x > for i = 1, . . . , n and ignoring the facets defined by these. 

A more efficient method (on some examples) is to find all facets and then 
remove the non-flippable irredundant facet normals by explicit checks. In our 
implementation this is done by checking if the inequality system with the in- 
equality in question inverted still has a positive solution. 

As mentioned in there is an algebraic test that helps us eliminate redun- 
dant inequalities of C^(I). Let a G K n be a coefficient vector of an inequality. 
If a indeed is irredundant and defines a facet with a relative interior point v 
then Corollary 12 . 141 tells us how to compute Q^(m v (I)). This marked reduced 
Grobner basis can be computed from Q^{I) as {^v(g)} g eg^(i) ^ we J us t know a 
and not necessarily v, see the next section. A necessary condition for a to be ir- 
redundant is that the computed set {^v(d)} g eg^(i) indeed is a marked Grobner 
basis i.e. all S-polynomials reduce to zero. This check even works for v outside 
the positive orthant. A quicker necessary condition that we can check is that 
every non-zero S-polynomial should have at least one of its terms in in^(I). For 
huge sets of inequalities the test works extremely well — 500 inequalities might 
reduce to 50 of which maybe 10 are irredundant. Our experience is that having 
this test as a preprocessing step can be much faster than solving the full linear 
programs with exact arithmetic. 

4.2 Local change 

Let Q^{I) be a known marked Grobner basis and let F be a flippable facet of 
C^(I). We let flip(C/^(I), F) denote the unique reduced Grobner basis different 
from G^(I) whose Grobner cone also has F as a facet. We will describe an 
algorithm for computing flip(<5^(/), F) given Q^,{I) and an inner normal vector 
a for F. For a marked Grobner basis Q and a polynomial g we let denote 
the normal form of g modulo Q and note that this form does not depend on the 
term order but only on Q. 

Algorithm 4.2 

Input: A marked reduced Grobner basis Q^{I) with -< being an unknown term 
order and an inner normal vector a of a flippable facet F of (I) . 
Output: g = flip(^(I),F). 

{ 

Let v be a positive vector in the relative interior of F; 
Compute G^(m v (I)) = {m v (g) : g £ £_<(/)}; 
Compute the marked basis £/^_ a (in„(I)) from Q^(m v (I)) 
using Buchberger's algorithm; 



18 



Q-={g-g GAI) :geg^_ a {m v (I))}; 

Mark the term m^_ a (g) in each element g — g&^i 1 ) in Q; 
Turn Q into a reduced basis; 

} 

The algorithm is a special case of the local change procedure for a single step in 
the Grobner walk [5]. See |HJ Proposition 3.2] for a new treatment and a proof. 
Here we will just add a few comments on our special case — the case where F 
is a facet and not a lower dimensional face: 

For any vector oj in the relative interior of F, in^(J) = in„(I) is homoge- 
neous with respect to the w-grading. Since F is (n — l)-dimensional, m v (I) is 
homogeneous with respect to all vectors inside span R (a) ± . All Grobner bases 
of in„(/) are homogeneous in the same way. Consequently, each of them must 
consist of polynomials of the form X^ =0 c s x^ a+sb ' where a G N n and i) 6 Z™ 
is parallel to a. The same is true for all polynomials appearing in any run of 
Buchberger's algorithm starting from one of these sets. A consequence is that 
in order to run Buchberger's algorithm we only need to decide if we are in the 
situation where x 7 -< x 7+a for 7 E N n or in the situation where x 7+Q -< x 7 for 
7 € N ra . Thus specifying a or —a as a term order suffices — no tie-breaker is 
needed. The initial ideal in v (I) can have at most two reduced Grobner bases. 
Both term orders are legal since in^(/) is homogeneous with respect to the 
strictly positive vector v. 

The Grobner basis £7^(in„(J)) can be read off from the marked Grobner 
basis Q~i,(I) by taking initial forms of the polynomials with respect to v, see 
Corollarv l2.14l Taking the initial form in v (g) of a polynomial g € G-<(I) without 
computing v is done as follows. By Corollary 12. 1H m^(in v (g)) = in^(g) and 
thus we already know one term of 'm v (g) since in^(g) is the marked term of 
g in Q^il). Since every oj in the relative interior of F will have t/^(in w (/)) = 
£/^(inu(/)) the remaining terms of in v (g) are exactly the terms in g with the 
same w-degree as in^(g) for all u in the relative interior of F and consequently 
for all oj in span IR (a)^. In other words a term of g is in in„(<7) if and only if its 
exponent vector minus the exponent of in^(g) is parallel to a. The term order 
-< does not have to be known for this step, nor does it have to be known in 
the computation of Q^_ a (in v (I)) or in any other subsequent step. The vector 
v also remains unknown in the entire process. 

4.3 Computing the search edge 

Let -< be the term order used for orienting the graph of the Grobner fan. In 
Algorithm 13.21 the search edge &q has to be computed given G^'(I) where 
-<' is some unspecified term order. According to Proposition 13. II the definition 
of search edges can be arbitrary. However, efficiently computing a search edge 
requires a good definition. Our search edges will always come from flippable 
facets. 

One strategy for locally computing the search edge ec%/(7) is to compute a 
unique representation of each flippable facet of the Grobner cone C^>(I) and 
then choose the smallest of these facets to be &c f (i) i n some lexicographic 
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order. This method requires all facets to be computed every time we check if 
"e u = (u,v)" in Algorithm EHl 

A better strategy is to draw a straight line from a point in the cone C^/(I) 
to the cone of the sink and choose the first facet intersecting this line as &c 
A point in the cone C^i(I) can be computed deterministically by linear pro- 
gramming. Two problems arise. The straight line might not intersect a unique 
facet and we may not know a point in the cone of the sink. Both problems 
can be solved using formal perturbation of the end points of the line. This was 
worked out in detail in [S] . Here we explain how it works for lexicographic term 
orders and with one end point perturbed. 

Lemma 4.3 Let I C R be an ideal and -< the lexicographic term order with 
x\ >- X2 >- ■ ■ ■ >~ x n . Define t £ = (e ,^ 1 , . . . ,e n ~ l ). There exists a 5 > such 
that in Ts (J) = in^(7) for all e 6 (0,5). 

Proof. This follows from Lemma 12.101 since -< and t £ agree on a finite set of 
polynomials for small e. □ 

Let a be a deterministically computed interior point of the cone of C^'(I) 
and assume for simplicity that a £ N n . For sufficiently small e > the line 
segment 

u(t) := (1 - t)a + tr £ with t € [0, 1] 

intersects a facet of C^(I) unless C^>(I) is the sink. 

Let {ai, . . . ,a m } be the set of potential inner facet normals read off from 
Q^i(I). We are only interested in the vectors where (cr, Oj) > and (t £ , oti) < 
0. Let ti denote the t-value for the intersection of the line segment and the 
hyperplane defined by Qj. Then 

= (cr, Oj) 

(a, - (r e ,ai) 

We wish to find i such that ti is smallest (for small e). 



ti 


< 




(6) 


(cr, a t ) 


< 


{o-,aj) 


(7) 


(cr, at) - {r £ ,ai) 


(cr,aj) - (T £ ,aj) 


(a, a t ) - (r £ ,ai) 


> 


(o-,aj) - (r £ ,aj) 


(8) 


(cr, a t ) 


(a,aj) 


(r E ,ai) 
(cr, Oi) 


< 


(T £ ,aj) 

(a,aj) 


(9) 


(t £ , (cr, aj)ai) 


< 


(T £ ,(a,ai)aj) <^ 


(10) 






x {<T,a i )a J 


(11) 



We see that for e sufficiently small '% < tj" does not depend on e. Furthermore, 
there cannot be any ties, unless ai and aj represent the same hyperplane. This 
gives an easy method for defining and computing ec< ,(/)■ We simply choose the 
facet defined by aj where U is smallest among {ti, . . . , t m } (for small e > 0). 
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4.4 Exploiting symmetry 

In this section we explain how to take advantage of symmetry to speed up 
computations. The symmetric group S n acts on polynomials and ideals of R 
by permuting variables and on R n by permuting coordinate entries. Let I C R 
be an ideal. We call a subgroup T < S n a symmetry group for / if tt(I) = I for 
all 7r € r. If we know a symmetry group for / we can enumerate the reduced 
Grobner bases of / up to symmetry. Let T be such a symmetry group for /. 

In our description all Grobner bases will be marked and reduced. Thereby 
each one will uniquely represent its initial ideal and Grobner cone. For a 
Grobner basis Q of / we use the notation Tg = {vr(^)} 7rg r for its orbit. 

The idea is to exploit the identity fup(7r(£?), tt(F)) = 7r(flip(C7, F)) for all 
7r £ r. In other words T is a group of automorphisms of the graph of the 
Grobner fan of /. The quotient graph is defined to be the graph whose vertices 
are the orbits of Grobner bases with two orbits Tg and Tg' being connected if 
there exists a facet F of the Grobner cone of Q such that flip(£?, F) € Tg*. The 
flip graph may have loops. 

The symmetry-exploiting algorithm enumerates the quotient graph by a 
breadth-first traversal. Orbits are represented by Grobner basis representatives. 
One question that arises is how to check if two Grobner bases Q and Q' represent 
the same orbit. A solution is to run through all elements ir E T and check if 
tt(G) equals Q' , or even better to make a similar check for the monomial initial 
ideals. Although this does not seem efficient, it is still much faster in practice 
than redoing symmetric Grobner basis and polyhedral computation as we have 
done in the usual reverse search or breadth-first enumeration without symmetry. 
It is not clear how to combine symmetry-exploiting with reverse search. 

5 Complexity 

In this section we will discuss the complexity of enumerating the maximal cones 
of the Grobner fan of an ideal / by reverse search. We will assume that / is 
homogeneous with respect to a positive vector. This guarantees that any facet 
of a full-dimensional Grobner cone is flippable. 

We identify the following important sub-algorithms: 

• Computation of the facet normals of the Grobner cone of a marked re- 
duced Grobner basis Q. We denote the time for this operation by Tf acets (G). 

• Computation of a search edge given a marked reduced Grobner basis Q 
as described in Subsection 14.31 We denote the time for this operation by 

Ishoot(<?)- 

• Conversion of a marked reduced Grobner basis Q\ into a marked reduced 
Grobner basis Q2 where the Grobner cones of Q\ and Q2 are assumed to 
share a facet. We denote the time for this operation by Tfn p (Gi,Q2)- 

For simplicity we will assume that the time used for solving any of the linear 
programs in our algorithms only depends on the dimensions of its matrix form. 
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We let Ti p (n,r) be the time needed to solve a linear programming problem with 
n variables and r constraints. Then Tf acets (C/) and T s h oot (G) can be expressed 
in terms of Ti p (n,r). 

The time Tf ace t s ({?)£ 0{Ti p (n,r)r) where r is the number of non-leading 
terms in Q. The reason is that each non- leading term in Q gives an inequality in 
the description of the Grdbner cone. Checking if the inequality defines a facet 
takes one linear program. In addition duplicates should be removed from the 
set of facet normals and further vectors should be eliminated until no parallel 
vectors exist. The time for this step is dominated by the time for solving LPs. 

The time T s h oot (G)£ 0(rn 2 +T\ p (n, r)) where r is the number of non-leading 
terms in Q. The first step in the algorithm is to deterministically find a relative 
interior point of the Grdbner cone. This is done in time Ti p (n,r). After this 
the smallest vector among the r defining vectors for the cone with respect to 
the ordering in Subsection 14. M needs to be found. Comparing two vectors takes 
0{n 2 ) operations in the worst case. These are operations in Q. In the above 
estimate we assume that each operation takes constant time. 

We have no good bound for the complexity of flipping. Now we count the 
number of times each of the three sub-algorithms are applied when enumerating 
the graph of the Grdbner fan of / using reverse search. 

• The facets of each Grdbner cone are computed exactly once in Algorithm 
13.21 (right after the Grdbner basis has been output). We remark that since 
we are only interested in facets with the correct orientation the number 
of LPs that really need to be solved is lower than the r in the discussion 
above. We will not take this into account in our analysis. 

• Checking if an edge is a search edge is done once for every edge. Every 
time we need to recompute a search edge and compare it to the edge. 
Hence the total number of times we need to compute a search edge is 
equal to the number of edges in the graph of the Grdbner fan of /. 

• When a vertex Q is processed by Algorithm 13.21 we must test for every 
ingoing edge if the edge is a search edge. To test this we first compute 
flip(C/,i ? ) where F is the facet of the cone corresponding to the edge 
in question. If the edge is equal to the search edge of Qip(Q,F) we do 
an enumeration of the subtree with root Rip(Q,F). If not, Rip(Q,F) is 
forgotten. Since all vertices are processed once and every edge is ingoing 
for exactly one vertex the number of times flip(£, F) needs to be computed 
is equal to the number of edges in the graph. We remark that the variant 
of the reverse search where the search path for the current vertex is not 
stored on the recursion stack would require twice as many computations 
of this kind. 

Let E be the edges and V be the vertices of the graph. The total time complexity 
of the enumeration of (V, E) is: 

0(J2 T facets(S) + Yl T shoot(£l)+ ^ T m P (Gl,G2}) 

S£V (ffi,&)e-B (£a,e 2 )e£ 
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Figure 4: The Grobner fan of the ideal in Example 16.11 intersected with the 
standard 2-simplex. The a-axis is on the right, the 6-axis on the left and the 
c-axis at the top. 

Substituting with the time needed for solving the LPs we get the following 
theorem: 

Theorem 5.1 Let (V,E) be the graph of the Grobner fan of I. The time 
complexity for computing this graph given a marked reduced Grobner basis is in 
the class of functions 

0(Y,Ti P (.n,r(g))r(g) + £ r{Q x )n 2 + T lp (n,0i) + £ Tm^, &)) 

GeV (Gi,G 2 )eE (G 2 ,Gi)eE 

where r(Q) is the number non-leading terms in the marked reduced Grobner 
basis Q . In particular, the first two terms are bounded by a polynomial in the 
size of the output. 

Corollary 5.2 If for a given class of ideals the time 7fli p (£/i, ^2) is bounded by 
a polynomial in the size of the binary encoding of Q\ and Qi then the enumer- 
ation of the reduced Grobner bases for an ideal in the class by reverse search is 
a polynomial time algorithm in the size of the output. 

6 Computational results and examples 

The algorithms presented in this paper were implemented in the software pack- 
age Gfan .12]. In this section we present examples of Grobner fans computed 
using this package. The first example comes with a picture and gives an idea 
of the kind of geometric shape a Grobner fan might have. 

Example 6.1 Q21 Example 3.9] Consider the ideal I = (a 5 - 1 + c 2 + b 3 , b 2 - 
1 + c + a 2 , c 3 - 1 + b 5 + a 6 ) C Q[a, b, c]. The Grobner fan of I has 360 full- 
dimensional cones and the Grobner region is M^q. The intersection of the fan 
with the standard simplex in 1R 3 is shown in Figure 0J 
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Example 


n 


h 


d 


D 


/-vector 


Det 3i3i 4 


12 


6 


3 


3 


(1,12,66,204,342,288,96) 


Det 3j 3 i5 


15 


7 


3 


3 


(1,45,585,3390,10710,19890,21750,12960,3240) 


Det3,4,4 


16 


7 


3 


5 


(1,?,?,?,?,?,?,?,?, 163032) 


Detsym 3 4 


10 


4 


3 


8 


(1,518,5412,20505,36024,29808,9395) 


Grass2,5 


10 


5 


2 


3 


(1,20,120,300,330,132) 


Cyclic 5 


5 





8 


15 


(1,?,?,?,?,55320) 


J 4 


4 


1 


3 


8 


(1,200,516,318) 



Figure 5: Statistics for the Grobner fans computed using Gfan. 

We now list some families of ideals used in our computations. The Grobner 
fans of these ideals have been computed for the parameters listed in the table 
of Figure El The ambient field is always Q. The columns of the table are to 
be interpreted as follows. In each row, the first column contains the name of 
the ideal (to be explained below). The second column lists n, the number of 
variables in the ideal. The third column lists h, the dimension of the lowest di- 
mensional Grobner cone Cq{I). Note that h is the dimension of the homogeneity 
space of the ideal which is the common subspace contained in every Grobner 
cone of the ideal. The quantity "d" is the lowest total degree of any reduced 
Grobner basis of the ideal and "D" is the highest. The f -vector of the Grobner 
fan is an ordered list of the number of /i-dimensional cones, h + 1-dimensional 
cones etc., up to the number of n-dimensional cones. 

Example 6.2 Let Dett jmjTl denote the ideal in the polynomial ring in ran vari- 
ables generated by the t x t minors of the matrix: 



X21 



X\2 
X22 



\ X m \ X m 2 



X\r, 
X2r, 



J 



Example 6.3 Let Grassy denote the ideal in the polynomial ring in Q) vari- 
ables generated by the relations on the d x d minors of a d x n matrix. 

Example 6.4 Let Detsym 4 n denote the ideal in the polynomial ring in ra ( TO 2 +1 ) 
variables generated by the t x t minors of a symmetric matrix of variables. For 
example, Detsym 3 4 is generated by the 3x3 minors of the following matrix: 



( a 


b 


c 


d\ 


b 


e 


f 


9 


c 


/ 


h 


i 


V d 


9 


i 


J ) 



Example 6.5 Let Cyclic 5 denote the ideal (a + b + c + d + e, ab + be + cd + de + 

ae, abc + bed + cde + dea + eab, abed + abce + abde + acde + bede, abede — 1) C 
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k[a, b, c, d, e\. In general, Cyclic n stands for the generalization of this polynomial 
system to n variables [Sj- These polynomial systems have become benchmarks 
for computer algebra packages and their lexicographic Grobner bases are noto- 
riously hard to compute. 

Example 6.6 Let K n denote the complete graph on n vertices and I n be the 
Stanley- Reisner ideal of this graph. The Stanley-Reisner ideal of a simpli- 



cial complex A is the ideal generated by all monomials x^Xi 



Xi„ 



such that 



{ii, . . . , i n } is not a face of A. Apply a generic linear change of coordinates to 
I n to obtain the ideal J n . The generators of J n typically have very complicated 
coefficients. For example, the first generator in our J4 was 



a 



3 I 4980248985 ^2 

T 34333RKfi4 



1980248985 2 „ 
343338664 u c 
359584197 r 2 
171669332 uc 
1194144Q14 ,2 j 
'"9656399925" " 
4583330213 j2 
3862559970 



I 2079196217 nh . 
257503998 uuu 
84523033581 3 
1373354656 L 
I 394500908221 . 
25750399800 ULU 
I 181743499 h j2 
^ 364392450 u 



+ 



86858380 
128751999 



b 2 ( 



11737327991 2 j 
51500799600" " 
47953955497 l j 
25750399800 uca 
429736138279 j2 
25750399800 LU 



2205648949 2 
42917333 uc 
16299027451 hrj 
38625599700 UU 
, 195491595943 2 j 
2985553600 c " 
I 8566043731 ]3 
~r 12875199900 " ' 



The initial ideals of J n are known as the generic initial ideals of I n . The reverse 
lexicographic generic initial ideals of an ideal have played an important role in 
commutative algebra and algebraic geometry while other generic initial ideals 
have not been explored too much. We computed the Grobner fan of J4. 

Extracting the f-vector from the full-dimensional Grobner cones produced 
in the enumeration process was the most time-consuming part of the compu- 
tation of these examples. In example Det3 3 4 this extraction was not possible 
to complete within reasonable time with the current software package. For this 
particular example the 163032 full-dimensional Grobner cones were computed 
up to the action of a symmetry group of order 576. The full-dimensional cones 
come in 289 orbits. The computation of the full dimensional cones up to sym- 
metry took 7 minutes on a 2.4 GHz Pentium processor. Using reverse search 
without symmetry the same computation would take approximately 14 hours. 
The f-vector extraction routine in Gfan only works for complete fans. This is 
why the f-vector for the Cyclic 5 example is not shown. 
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